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The joint action of a synthetic gauge potential and of atomic contact repulsion in a one-dimensional 
alkaline-earth(-like) fermionic gas with nuclear spin I leads to the existence of a hierarchy of frac¬ 
tional insulating and conducting states with intriguing properties. We unveil the existence and the 
features of those phases by means of both analytical bosonization techniques and numerical meth¬ 
ods based on the density-matrix renormalization group algorithm. In particular, we show that the 
gapless phases can support helical modes, whereas the gapped states, which appear under certain 
conditions, are characterised both by density and magnetic order. Several distinct features emerge 
solely for spin / larger than 1/2, thus making their study with cold-atoms unique. We will hnally 
argue that these states are related to the properties of an unconventional fractional quantum Hall 
effect in the thin-torus limit. The properties of this hierarchy of states can be experimentally studied 
in state-of-the-art cold-atom laboratories. 


The simultaneous presence of particle-particle inter¬ 
actions and of (non-)Abelian gauge potentials, such as 
magnetic fields or spin-orbit coupling (SOC), is respon¬ 
sible for several spectacular phenomena, the fractional 
quantum Hall effect (QHE) being only the most known 
example [1]. Since several years this interplay is under 
close scrutiny both because of its perspective role in 
the realisation of robust quantum information protocols 
and because of the general interest in topological states 
of matter [5]. In the presence of SOC [2 0], interac¬ 
tions can further drive the system into fractional quan¬ 
tum spin Hall states mi where edge-currents are spin 
polarised, or other exotic phases, characterised by un¬ 
usual spin textures [5], fractional conducting modes [5] 
or parafermions nni. 

Up to now, the attention has almost entirely focused 
on spin-1/2 (electronic) liquids, as most appropriate for 
the description of condensed matter systems. However, 
the recent progresses in the manipulation and control 
of cold atomic gases m [12] have brought to high rel¬ 
evance the study of systems of interacting fermions with 
a large (and tunable) spin. The investigation of alkaline- 
earth(-like) atoms such as Ytterbium [TSHTT] or Stron¬ 
tium usms], which are characterised by a nuclear spin I 
larger than 1/2, is opening the path to the exploration of 
phenomena, e.g. the properties of SU(A/’) models [20ll^ . 
which are not accessible with solid-state systems. The 
phase diagram of related multi-component Heisenberg 
or Hubbard-like models has been investigated in several 
works [221 - I?7] : yet, very little is known concerning the 
effect of a gauge potential on an interacting system of 
particles with a large spin (see, however, Ref. [28]). 

Synthetic gauge potentials in cold atomic systems can 
be induced via properly tailored laser pulses [35] . The im¬ 
plementation of these schemes has already led to the re¬ 
alisation of light-induced magnetic fields [30] , to Rashba 
SOC [3T] , as well as to lattice models with non-zero Chern 


numbers [32H37] and ladders with synthetic gauge poten¬ 
tials [35] • Similar approaches are suitable for application 
also to multi-component gases iMiiin], as shown in two 
recent spectacular experiments mmi]. This fecund ex¬ 
perimental activity, together with the rich scenario al¬ 
ready explored for spin-1/2 systems, motivates the inves¬ 
tigation of interacting systems with large spin coupled to 
gauge potentials. 

Do large I alkaline-earth(-hke) atoms lead to mere ex¬ 
tensions of what is already known for electronic liquids? 
The answer is no. We will provide examples confirm¬ 
ing that these setups allow for the exploration of novel 
regimes that can naturally be achieved only for / > 1/2 
or through the versatility of this new playground. This 
puts cold-atom experiments in an excellent position for 
the investigation of novel intriguing many-body effects 
unattainable in conventional condensed matter setups. 

In this article we consider a one-dimensional fermionic 
gas with nuclear spin / > 1/2 and investigate the joint 
effect of interactions and of a synthetic gauge potential, 
which is equivalent to a Rashba SOC and an external 
magnetic field. Provided that the states with highest and 
lowest spin are directly coupled with a multi-photon tran¬ 
sition im SO], a full hierarchy of magnetic crystalline 
states appears at fractional fillings: 

TTTh p 

^ = 1- tor I 1 ^ = “! p, g G N+ and co-prime (1) 

fcso (2/ + 1) q 

where n is the atomic density and fcgo is the typical 
momentum of the SOC (to be defined in the following). 
Combining analytical and numerical methods, we show 
that these insulating phases exhibit charge and spin or¬ 
dering and in some cases are connected to an unconven¬ 
tional fractional quantum Hall states in the thin-torus 
limit [43] , while intrinsically resembling an edge of a spin 
Hall system. The stabilisation of these gapped phases 
with g > 1 requires some form of atom-atom interaction. 
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Whereas some of the phases can be realised in the pres¬ 
ence of a simple contact repulsion, it is in general true 
the higher the g, the longer the range of the necessary 
interaction. These phases change dramatically once the 
mentioned multi-photon coupling is switched off: simi¬ 
larly to the spin-1/2 case [5], fractional helical liquids 
appear. 

Our analysis, which includes also the effect of a 
trapping potential, confirms that the findings of this 
work can be observed with state-of-the-art experimen¬ 
tal techniques. Following the reasoning put forward in 
Ref. [3^130], we can conclude that this setup might serve, 
especially in the limit of large /, as a quantum simulator 
of two-dimensional exotic interacting phases of matter in 
the strongly anisotropic limit. 

Model. We consider a one-dimensional model that 
describes an optical lattice loaded with a gas of fermionic 
alkaline-earth(-like) atoms whose ground state is char¬ 
acterised by a spin-/ nuclear manifold; for a sketch, see 
Fig .[^a). Disregarding for the moment the harmonic con¬ 
finement, the Hamiltonian reads [50] 

I 

Tdo = + H.C.^ + "Hint ; (2) 

j m=-I 

where are fermionic operators annihilating (creat¬ 
ing) an atom at site j with nuclear spin m and t is 
the hopping amplitude. Hint describes an SU(2/ -|- 1) 
invariant interaction which is usually of contact kind: 
^ ■ As new proposals make the en¬ 

gineering of longer-range interactions in cold gases more 
practicable |33|) we also discuss as an example the effect 
of a nearest-neighbour potential: V m m' 

A Raman coupling endowed with a running phase con¬ 
nects states which differ for one nuclear magnetic quan¬ 
tum Am = ±1: 

/-I 

= E E + H.c.) ; (3) 

3 rn=-I 

here, Vim = with D the Raman-coupling strength 

and gm = \JI{I -\-1) — m{m + 1); see Ref. |30| for a 
derivation of 

The unitary transformation U defined by = 

maps Hamiltonian Hq -1- Hi to a spin-/ 
fermionic model in the presence of Rashba SOC and 
of a magnetic field H with perpendicular quantization 
axis (see Supplementary Material). The choice to de¬ 
note the phase factor in Eq. (§ with /so becomes then 
clear upon inspection of the kinetic term: U Hq W = 
2t ^^k.m ^Os(/ 277iA:sOA Hint- 
Alternatively, Hq + Hi can be interpreted as a spin¬ 
less fermionic model with one (finite) synthetic dimension 
coupled to a synthetic Abelian gauge field m- Multi¬ 
photon transitions can implement an additional coupling 





FIG. 1. One-dimensional gas of^^^Yb atoms, (a) Sketch of a 
one-dimensional gas of atoms with / = 5/2. (b) Energy 

bands of Hamiltonian Ho for Hint = 0 and I = 5/2; an energy 
shift is inserted for representation clarity. When Q is turned 
on, fermions with momentum difference Afc = 2/so and spin 
difference Am = ±1 (solid arrows) or Am = ±2/ + 1 (dashed 
arrows) get coupled through Hi and Ho, respectively: if the 
condition kso = kp is met, the system develops a full gap, cor¬ 
responding to = 1. (c) For 1 = 1/2 the condition kso = kp 
is not enough because Ho = 0 (upper panel): when kso = 7r/2 
the identification of momenta modulo 27r allows for the cre¬ 
ation of a gap (lower panel), (d) When Hint yf 0, the system 
can develop a gap for lower fillings v = l/q via higher-order 
scattering terms. As an example, the picture highlights three 
intermediate processes: their sequence (top to bottom) origi¬ 
nates a third-order process which couples two Fermi surfaces 
with Am = ±1 for q = 3 and kp = kso/3. 


between the extremal nuclear spins m = ±/, described by 
a contribution Ho = -|-H.c.^ to 

the Hamiltonian |45| . In the following we consider Yt' = Yl 
because, especially for small /, multi-photon transitions 
Yl' can be tuned to the order of magnitude of the Raman- 
coupling Ylm. The presence or absence of Ho corresponds 
to different boundary conditions along the synthetic di¬ 
mension, a feature that is relevant only for / > 1/2 be¬ 
cause, for / = 1/2, Hamiltonian Ho connects states that 
are already linked hy Hi. The possibility to tune (and 
eventually switch off) Ho is unique to this cold-atom im¬ 
plementation. 

Finally, the presence of a trapping potential is taken 
into account through an additional term in the Hamil¬ 
tonian of the form 7/tr = m '^3 mH,m- The model 
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FIG. 2. Charge and magnetic order of the magnetic crystals. Charge (hi) and magnetic order (M“) of the fractional phases at 
I and I (from left to right) for 1 = 1, fcso = 7r/3 (a) and I = 1/2, kso = 7r/2 (b) as obtained from DMRG simnlations 
of a system of length L = 96 with Q/t = 1. For the interaction parameters see the panels. Since the system is a crystal with 
small boundary effects, we only plot the central part of the system for a better readability, (c): Sketch of the density (red balls 
with radins related to {nj}) and magnetic properties (blue arrows representing the vector (M“)) of the insulators at v = 1, | 
and I (from left to right) for 1 = 1. 


defined by the Hamiltonian Hq + Hi + ^.2 + Htr will 
be studied in its various limits by means of bosonization 
techniques |46) and density-matrix renormalization group 
(DMRG) simulations [47l |48j , which are specially suited 
for interacting systems in one dimension. 

Magnetic crystals. As a first step we show that 
for the fillings given in Eq. Q the system described by 
Hamiltonian Hq +711+ ^2 and with / > 1 is insulating. 
Note that there is no inconsistency between the fermionic 
statistics and the study of integer I, as the considered 
states can be selected from a larger half-integer mani¬ 
fold [iniUn]- The peculiar case I = 1/2 is also considered. 

Gapped phases appear in the absence of interactions 
only for the fractions in Eq. ([^ with q = 1. It is in¬ 
structive to briefly discuss this case as it helps under¬ 
standing the interacting one. By writing Hamiltonian 
Tio + Hi + H 2 in momentum-space representation and 
defining a Fermi energy and a Fermi momentum kp rel¬ 
ative to the case with 11 = 0, it is immediate to see that 
whenever fcso = kp, a. fermion can be spin-flipped and 
scattered from one edge of the Fermi sea to the other 
one. As any edge of the Fermi sea gets coupled, the sys¬ 
tem develops a full gap (see Fig. [^b)). In this case, the 
density n is (2/ -|- l)kpjTr, so that p = 1. Because with 
higher-order processes Hi and H 2 can connect spin states 
with |Ato| > 1, a gap opens whenever kso = kp/p, 
p G [1, 2 ... 2/]: these phases correspond to p = p. 

The case I = 1/2 is peculiar since H 2 = 0 and only 
two edge-modes gap out [33]. Exploiting the presence of 
a lattice, it is however possible to invoke the equivalence 
of momenta upon shifts of 27r and hence to derive the 
additional condition kp + 2kso = —kp + 27r. Together 
with kso = kp (in this case kso = it ensures that 

the system enters a gapped phase aip = 1 (see Fig. Sc)). 


Contrary to the previous case, the resulting gapped phase 
is a lattice effect, without a proper continuum limit. 

The existence of gapped phases for 9 > 1 and odd can 
be predicted using the bosonization technique [461IM1I5T] . 
which is the natural analytical tool to take into account 
interactions in one dimension (see the Supplementary 
Material for a discussion of the case of even q, together 
with more details about the bosonization approach). By 
linearizing the non-interacting spectrum of Ho around 
the Fermi energy, we can write Cm{x) = -b 

where (V'-.m) is the right (left)- 

moving operator of the m-th nuclear spin state. Central 
to the theory is its expression in terms of the bosonic 
fields 4>r,mix), r = ±, as 'ifr,7n{x) ^ exp [-jr 0 r,m(a;)], 
with [(i)r,m{x),^r' ,m'{x')] = -i r TT Sgn(x' - x)5ry Sm,m'- 
In the presence of density-density interaction terms, such 
as Hint, the full Hamiltonian Ho can be cast into the 
quadratic form 

■^0 = X! / d.^idx+mix) dJm{x))Mm,m' 

m,m' ^ ' 

(4) 

with Mm,m' = vp/{2TT)6m,m' +Um,m'', vp is the Fermi 
velocity, and 4’r,m{x) = (pm{x) — r6ra(x)',lAm,m' describes 
the scattering processes induced by Hint only. 

Insulating phases are determined by the joint action 
of Hi + H 2 and Hint, representing the interplay of the 
gauge potential and interactions. Indeed, Hi originates 
21 processes of the form 

,"1+1 ,"i+l) J G N ; 

(5) 

with m = — 1 , which are relevant only when 

momentum is conserved, namely kp = kso/Q with q = 
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2n + l. When n = 1, for instance, term ([^ reproduces the 
low-energy physics of the processes displayed in Fig.[^d). 
7^2 generates one more relevant process, expressed in the 
form of Eq. ([^ by replacing m with / and to-|- 1 with —I. 
Once Eq. (§ is written in terms of the bosonic fields, 21 -b 
1 Sine-Gordon Hamiltonians with commuting arguments 
are obtained. Following Ref. [sniEi], when these terms 
are relevant in a renormalization group sense, they lead 
to the formation of 2 / -b 1 mass gaps which make Tio + 
iii+ii 2 fully gapped. 

The request for relevance in the renormalization group 
sense points out an important fact: given a specific fill¬ 
ing, a generic interaction does not necessarily stabilize a 
gapped phase. Through the explicit study of the 7=1/2 
case, for which precise mappings between the microscopic 
interacting model and the bosonization parameters are 
known, in the Supplementary Material we argue that the 
higher the value of g, the longer the range of the inter¬ 
actions required to make the 27 -b 1 Sine-Gordon terms 
relevant. 

In order to fully characterise the properties of this hi¬ 
erarchy of gapped phases, we perform numerical simula¬ 
tions with DMRG [13 EH] that in one dimension provide 
essentially exact results and allow us to explore the cases 
of even q and p > 1 , which are not easily accessible with 
bosonization (see the Supplementary Material for details 
on the numerical simulations). We particularly focus on 
7 = 1/2 and 7=1 (for these cases gm = 9 and can be 
set to 1 ) and use the latter case to highlight the general 
new features emergiM for 7 > 1/2. 

The results in Fig. ^consider the fillings *^= 5,5 and | 
and display the density profile (fij) = and the 

magnetization (A7“) = , where 

S'" is a (27 -b l)-dimensional representation of the SU(2) 
spin operator {a = x, y, z). Fig. [^clearly shows that the 
incompressible phases are characterised by both charge 
and magnetic order. 

Let us begin with the charge order, considering for ex¬ 
ample the case 7 = 1/2, v = ^ \ here, a nearest-neighbour 
interaction stabilises a density wave with one-particle ev¬ 
ery three sites. Note that in the usual SU(2) Hubbard 
model (H = 0 , no gauge potential) a nearest-neighbour 
interaction only stabilises a density wave with one parti¬ 
cle every two sites. Conversely, for U,V = 0 (no interac¬ 
tions) no gapped phases exist at this filling. The simulta¬ 
neous presence of interactions and of the (commensurate) 
SOC proves to be crucial in all the considered cases both 
for the opening of the gap and for crafting the properties 
of the resulting insulator. 

Concerning magnetic order, in all the considered cases 
(Mj) = 0 , which implies that in every site the ±m 
nuclear-spin components have the same occupation num¬ 
bers. Magnetization lays in the x — y plane and winds as 
a function of position in peculiar ways which depend on 
the filling (see the sketches in Fig. for 7 = 1). For 
7 = 1/2, the Hamiltonian is real and thus (A7j) = 0: the 


magnetization only develops along the x axis. 

The density and magnetic ordering found in DMRG 
simulations becomes transparent by considering the limit 
Vl/t ^ 1, H/f7 1 where a simple model catches its 
salient features. Hamiltonian ?7i -b 772 is invariant un¬ 
der the Z 2/+1 group related to the on-site transformation 
m I—>■ TO -b 1 (here, 7 -b 1 = —7), and is diagonalized by 
27 -b 1 local eigenmodes of the form: 

= ^ ( 6 ) 

^ m 


with A = 0 ... 27. As 77o is SU(27 -b 1) invariant, it is un¬ 
affected by this basis change. However, the on-site energy 
of each transformed spin state depends both on the site 


j and on kso'- £j,\ = 2Hcos 


2ttX 

2/+1 


2A:soj 


. When the 


atomic density n and fcso are commensurate, the periodic 
energy landscape £j^\ induces charge-density waves, with 
atoms magnetised accordingly. As shown in Fig. this 
pattern ultimately determines the charge and spin order 
even at H/t ~ 1. For 7 = 1 and fcso = 7r/3, the system 
dimerizes (this is because £j^\ has two degenerate lowest- 
energy states) and for Jz = 1 and ^ dimers arrangements 


are quite apparent. For 7 = 1/2, the modes dj^\ are the 
eigenstates of MJ, and this explains why magnetization 
develops only along x. This simple model provides a phys¬ 
ical intuition of the fact that 77i -b 7/2 induces magnetic 
crystals at the fillings in Eq. Q, thus complementing the 


a 


b 



Jmod3=l j mod 3 = 2 j mod 3 = 0 jmod2=\ j mod 2 = 0 


c 




site(j) site[j) site(j) 


FIG. 3. Magnetic crystals. Energy-spin structure for: (a) 
7 = 1 and kso = 7r/3, (b) 7 = 1/2 and kso = 7r/2. Density 
plots (uj.a) = {dl xdj^x) at i/ = |, I and | (from left to 
right) for 7 = 1 (c) and 7 = 1/2 (d) obtained through DMRG 
simulations (see the caption of Fig. [^ for the parameters). 














































5 



FIG. 4. Harmonic confinement effects. Density profile of the 
system / = 1 in the presence of a harmonic confinement with 
w = 10“^. The parameters of the simulation, done for 16 
fermions, are U/t = 10, Q/t = 1 and kso = vr/3. In the center 
of the system the typical charge and magnetic order of the 
insulator with filling = | appear. 

information provided by bosonization for odd q and by 
DMRG for the other considered cases. 

The gapped phases in Fig. at :/ = ^ and | are sta¬ 
bilized by the solely on-site repulsion, and thus poten¬ 
tially experimentally accessible in current experimental 
setting. It is therefore important to check whether the 
presence of the trapping potential will modify our find¬ 
ings. This is explored in Fig. where we show that the 
presence of a harmonic conhnement, introduced through 
T^tr with Wj = tc (j — L/2 — 1/2)^ does not hinder the 
possibility of observing the peculiar properties of the 
magnetic crystals discussed so far; in a Thomas-Fermi 
spirit, they form in definite regions of the trap. It is nat¬ 
ural to envision that the peculiar density-spin patterns 
which characterise these gapped phases could be unam¬ 
biguously revealed through spin-resolved single-site ad¬ 
dressing [52l |53]. Alternatively, Bragg scattering could 
get access to the structure factor of the gas, which con¬ 
tains information about the density order and, if spin- 
resolved, even about the magnetic one |54j . 

Helical liquids. According to the previous discus¬ 
sion, for / > I gapped phases at the fillings in Eq. Q 
arise in the interacting system only when 'H 2 ^ 0 and, for 
7 = 1/2, when kso = (lattice effect). More generally, 
explicit inspection has shown that even when 7^2 = 0 the 
condition kso = (27 -I- 1) for 7 > 1 can open a full gap 

through a high-order process which exploits the presence 
of a lattice (i.e. of momenta identification upon shifts of 
27r). When these conditions are not met, and thus the full 
gap does not develop, we will show that for the fillings 
in Eq. 0 the system described by 77 0 -I- "Hi is a helical 
liquid 0 . 

An intuition for this phenomenology is again provided 
by the non-interacting case 77int = 0 and q = 1. Here, 
the two Fermi edges at m = ±7 and k = ikp, respec¬ 


tively, remain unperturbed (see Fig. [^b)) and represent 
two gapless helical fermionic modes, which are the lowest- 
energy excitations of the system. Bosonization techniques 
applied to the case 77int 7 ^ 0 for q > 1 and odd clearly 
pinpoint the helical nature of the low-energy spectrum, 
although the gapless modes for g > 1 are linear com¬ 
binations of the original modes. Additionally, their con¬ 
ductance is fractional [9]: these two properties define a 
fractional helical liquid. Similarly to the magnetic crys¬ 
tals, bosonization reveals the existence of requirements 
on the range and intensity of the atom-atom repulsion. 

Numerically we cannot fully access the helical nature of 
the first excitations; we rather diagnose the existence of 
a helical ground state through two observable quantities. 
First, as a consequence of the existence of two gapless 
modes, the low-energy spectrum is described by a con¬ 
formal field theory in 1 - 1-1 dimensions with central charge 
c = 1. Second, the system is characterised by a current 
pattern where = —z7c]_^Cj+i,m+H.c., which 

is different from zero and features a flow direction related 
to the sign of m. 

Our numerical results for the case 7=1, kso = tt/S 
and 1 / = I are shown in Fig. (for other fillings and 
values of 7 see the Supplementary Material). For ground 
states of theories with a low-energy conformal limit, the 
Calabrese-Cardy formula [5S] predicts that the entangle¬ 
ment entropy S{£) = —Tr[p^logp^], obtained through 
a bipartition of the system into two blocks of length 
£ and L — £ {pi is the reduced density matrix of the 
block of size £), should be proportional to c: S{£) = 
A + c/61og[2L/7rsin(77r/L)] (for open boundary condi¬ 
tions). Figure shows S(£) as computed through DMRG 
simulations for L = 192 and a fit with the previous for¬ 
mula leaving A and c as free parameters. The result 
clearly indicates c = 1. The current profile (77/^) dis¬ 
played in Fig.j^shows a hnite saturated value for j ~ L/2 
and m = ±1 (note that they flow in different direc¬ 
tions); for m = 0 the absence of a current is imposed 
by symmetry reasons. The oscillations are a boundary 
and finite-size effect and decay to zero for L ^ 00 (see 
Supplementary Material). Note that in this case, even 
if kso = 7r/(27 -|- 1 ), the system is gapless and thus 
amenable to the study of helical properties. 

The simultaneous presence of two gapless modes (c = 
1) and of helical current patterns, indicating the presence 
of a helical state, results from the interplay between in¬ 
teractions and the gauge potential. In the non-interacting 
case, at this filling the system has three gapless edge 
modes, which reflect in a central charge c = 3, as shown 
in Figure]^ On the other hand, if we tune fl = 0 (and 
thus kso disappears from the Hamiltonian), the model 
is equivalent to a SU(3) Hubbard model, which cannot 
develop helical currents, since the Hamiltonian is real. 
Even the milder request D 77 0, kso = 0 results in a state 
without helical currents for the same reason. 

Finally, we would like to comment on the role of 772 
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in the / > 1 case: the presence of the multi-photon tran¬ 
sition leads to the absence of gapless fractional helical 
modes and vice versa. Upon interpreting our setup as 
a spinless system with one synthetic dimension [351 SO] 
and boundary conditions depending on 'H 2 , we have a 
phenomenology reminiscent of that of topological mod¬ 
els with edge modes, which gap out once the boundary 
conditions are changed. In the non-interacting limit, this 
analogy becomes exact for fillings with g = 1 in the 
/ —> 00 limit, where the physics of the integer QHE is 
recovered. It is remarkable that this analogy holds also 
for interacting systems (g > 1) and for small I. This 
rich phenomenology is lost for 1=1/2 (or for spin-1/2 
quantum liquids, such as electron degenerate gases): in a 
two-leg geometry, there is no difference between a cylin¬ 
der and a stripe. 

Conclusions. Alkaline-earth(-like) gases provide 
access to the physics of (quasi-)two-dimensional mod¬ 
els through the mapping of the nuclear-spin states to 
a (finite) synthetic dimension |39j . In the system con¬ 
sidered in this article, the synthetic lattice is pierced 
by a magnetic gauge potential with flux per plaquette 
<i>/$o = 2A:so is the quantum of flux) [40], whose 
properties are thus related to the original SOC potential. 

Interestingly, for / > 1, the gapped phases that we have 
studied are amenable to an intriguing interpretation by 
introducing an extra synthetic-dimension. Indeed, when 
7/2 0 the effective geometry of the system is that of 

a narrow cylinder since the synthetic-dimension length, 
2/ -|- 1, is much smaller than the real-space one. In this 
limit, usually called the thin-torus limit (TTL) [33115fij . 
the fractional quantum Hall states are (continuously) 
turned into charge-density waves, with features which 
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FIG. 5. Helical liquids. DMRG simulations of the entangle¬ 
ment entropy of the state for 7 = 1, fcso = 7r/3, i/ = |, 
Q/t = 1 for U/t —>■ oo, L = 192 (a) and U/t = 0, L = 300 
(b) . Thin red lines are fits with the Galabrese-Cardy formula 
which yield c = 1.0 ± 0.2 (a) and c = 2.94 ± 0.06 (b). (c) 
Helical currents for fhe case Q = t. For clarity, only 

the first half of the chain is plotted. 


coincide with those presented in this article. First, the 
filling v coincides with the ratio between the number of 
particles and the number of magnetic fluxes piercing the 
synthetic lattice, iV/iV$, which is the well-known condi¬ 
tion for observing the Laughlin series in QHE. Second, 
the density waves found in the TTL for v = 1/q dis¬ 
play one particle every g sites, which is what we show 
in Fig. [^ Additionally, the helical properties which ap¬ 
pear for 772 = 0 can be interpreted as precursors of the 
edge-modes of the QHE [30] . 

It is remarkable to observe the presence of these fea¬ 
tures in our system, since the SU(2/ -|- l)-invariant inter¬ 
action is, in the quasi-two-dimensional picture, strongly 
anisotropic, short-ranged along the chain and infinite- 
ranged in the synthetic direction, and does not resemble 
the features neither of the Coulomb nor of the contact 
repulsion which are usually considered in the QHE the¬ 
ory. 

Spin-orbit coupled alkaline-earth(-like) gases are thus 
a natural quantum simulator of the physics of the frac¬ 
tional QHE in an array of quantum wires [51 [Sj. More 
precisely, because of the unusual properties of the inter¬ 
action, we are dealing here with an unconventional form 
of QHE, and it is an exciting perspective to investigate up 
to which point it shares features with the standard QHE. 
More generally, large-/ Fermi gases provide a valuable 
experimental toolbox for the study of two-dimensional 
exotic phases of matter through coupled arrays of one¬ 
dimensional systems [51 El EZHSl- Additionally, the 
Raman coupling which we have described may realize 
a Za-invariant model: a symmetry-class which supports 
topological phases with edge parafermions mm- 

Concluding, in this article we have provided an analyt¬ 
ical and numerical study of spin-orbit coupled alkaline- 
earth(-like) gases and demonstrated the existence of a 
full hierarchy of gapless and gapped phases with exotic 
properties. In particular, we have identified and char¬ 
acterised the phases which can be experimentally re¬ 
alised with state-of-the-art cold-atom experiments, rang¬ 
ing from magnetic crystals to helical phases [5], identi¬ 
fying the effect of the simultaneous presence of interac¬ 
tions and gauge potentials. Our study also suggests that 
alkaline-earth(-like) atoms are a promising tool for bring¬ 
ing ultra-cold atomic gases into the quantum Hall regime, 
a long-standing and yet to be achieved goal, through the 
access of its thin-torus limit. The unprecedented versa¬ 
tility of these setups motivates further speculations and 
research: we leave for example as an interesting open per¬ 
spective the extension of this study to multicomponent 
bosonic systems [53H35] . 
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Supplementary Material 

Magnetic crystals and helical liquids in 
alkaline-earth fermionic gases 


Relation with Rashba SOC. Let us first explicitly 
display that the unitary transformation lA defined in the 
Results leads to a Hamiltonian which is formally equiva¬ 
lent to a Rashba SOC model: 

U [Ho + Hi] = -t ^ + H.c.) + 

I-l 

“bHint “b ^ “b H.C.^ . 

j m=-I 


We have already commented on the fact that the term 
proportional to t is the lattice version of a Rashba SOC. 
Concerning interactions, it is natural to assume that Hint 
depends only on the density operators which are 

left unchanged by d; thus: U Hint = Hint- Finally, the 
term proportional to U is the result of a magnetic field 
applied along the x direction, perpendicular thus to the 
quantization axis of the SOC. In an open chain, such as 
in our numerical simulations, these avoids the possibility 
of gauging away the phase 2 fcso- 

Bosonization. In the following we briefly discuss the 
existence of the magnetic crystals within the bosoniza¬ 
tion framework. The discussion follows the guidelines set 
in Ref. [5J [SOI IS] • To this aim we formally rewrite the 
fermionic operator as 

c™{x) = 

n—0 


in order to take into account the non-linearities of the 
free spectrum which play a non trivial role when inter¬ 
actions are strong enough; are unknown coefficients. 
Applying (SI I to a single term c^^(^x)cm+i{x) of 

Hi or to H 2 we get several contributions, of which those 
in the form: 


(xL+i)”'; 

(S2) 

are particularly interesting. Here, we have introduced the 
notation: Xrn = xn = —I ,...,/ (with / -|- I = 

—I). The terms (S2| are the only ones which conserve 
momentum independently on the value of kgo^ provided 
that 


fej’= fcso/?! with <7 = n-b n'-f 1. (S3) 

When n = n' (and thus g = 2n -|- 1 is an odd integer), 
the bosonized version of such operators is: 


(see the Results section for the definition of the (^±,m(a;)). 
Note that [0^,0^.'] = 0 for m,m' = —I,..., I. When 
the operators Om {rn = —/,...,/) are relevant in the 
renormalization group sense (here we assume to be in 
the condition for which this is true), they are pinned 
to minimize the quantity cos[(!Im] and in a semi-classic 
approach, i.e. cosOm^rn+i « 1 — gapless 

Luttinger Hamiltonian Ho acquires 2/ -|-1 mass gaps (i.e. 
it becomes fully gapped). 

The existence of fractional phases with q even cannot 
be straightforwardly explained using bosonization. Let 
us consider as an example the case q = 2. For each term 


(S2) we can choose n = 0 and n' = 1 but also n = I and 


n' = 0 , so that 2 ( 2 / -|- 1 ) sine-Gordon terms with non 
commuting arguments and with the same scaling dimen¬ 
sion appear. In this case a semi-classic approach cannot 
be used and the solution of the resulting Hamiltonian is 
an interesting question. 

As a last remark, we mention that the bosonization 
framework can be used to show the existence of fractional 
insulating phases for p > 1 and q odd, by introducing a 
fictitious coupling between different nuclear spin states 
of the form m- 

i-p 


H1-HH2 

j Lm=-/ 






H.c. 


(S5) 


which models higher-order couplings between spin states 
with Am 7 b ± 1 . Bosonization of this term yields results 
which are completely equivalent to the previous ones, 
apart from the fillings at which gapped phases appear. 

When we consider atoms with two nuclear spin states, 
namely I = 1/2, the Hamiltonian H 2 = 0. Gapped phases 


at fillings v = 1/q with odd q are stabilised by (S 2 ) and 
by 


^2i[kF( — n—n' — l) — kso]x 


{Xr. 


’i’— + ,ra+l iXm+l) 


(S 6 ) 

which also appear in the previous expansion. With m = 
— 1 / 2 , they imply kso = 7J'/2 and kp = TT/(2q) re¬ 
spectively. As already discussed in the main text the 


term (S 6 ) disappears in a continuum limit. If we intro¬ 


duce the charge (c) and spin (s) bosonic fields (f>r,a- = 
1 /a/ 2 [i^c - r9c -b a{ips — ^^s)] the Hamiltoninan Hq can 
be written as 


X—c,s 


dx ( Kx{dJ\f -b 


(S7) 


cos 


+ '/-,m + 4>+,m+l + </-.m-|-l) + 


4^+,m “b 0 —,m+l 


= COS Or, 


(S4) 


K\ is the usual Luttinger parameter which takes into 
account the strength of the interaction, u\ = vpjKx is 
the renormalized Fermi velocity. On the other hand the 
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terms (S2) and 

?ii oc ^ J dx cos[V2{q<fc + V^s)] 


(S 8 ) 


77 =± • 


couple the charge and the spin degrees of freedom. A 
renormalization-group calculation shows that they are 
relevant when < 3/q^, assuming a SU(2) invari¬ 
ant interaction, i.e. Kg = 1 [H]- Physically, this means 
that an on-site repulsive interaction, for which 1/2 < 
cannot stabilize a gapped phase for q > 3. 
Longer-range interactions are thus necessary (for nearest 
neighbor repulsion, for instance, it is possible to achieve 
Kc < [7D]), and in general we expect that, even for 

/ > 1, the higher the value of g, the longer the range of 
the interactions required to open the gap. 

Numerical simulations. DMRG is an algorithm 
which performs a search of the ground-state of a Hamil¬ 
tonian in the space of matrix-product states (MPS), a 
class of states with finite correlations characterised by 
the so-called bond link, D [171 US]. In the limit D = 1 
MPS are product states, whereas for larger values of D 
more quantum correlations can be described. 

We consider chains with open boundaries and length 
comprised between L = 96 and L = 192; setting D = 150 
we are able to describe the correlations in the states with 
sufficient accuracy. With these parameters the effect of 
boundaries is irrelevant and the errors on the observables 
are negligible on the scale of the symbols used in the fig¬ 
ures. In the simulations of the gapped crystals conver¬ 
gence is helped by the quantum numbers related to the 
conservation of each magnetization moreover we 

find that it is important to alternate the infinite-size ver¬ 
sion of the algorithm with the finite-size one. 


Supplementary Material for the case 7 = 1/2 

Numerical results for kgo 7 r/ 2 . Let us numer¬ 
ically show that fcso = 7r/2 is a necessary condition 
to obtain crystalline phases for I = 1/2. We con¬ 
sider as an example v = 1/2 and the interaction 77 = 
=i/ 2 ^j,m=-i /2 in the limit U /t — >■ oo, that sta¬ 
bilises a gapped phase for kso = 7 r/ 2 . Consistently, the 
entanglement entropy of the ground state displays a clear 

a) . If kso is rather tuned 


area-law behaviour, see Fig. SI 


to 7 r/ 3 , under the same conditions the entanglement en¬ 
tropy shows a non-area-law behaviour, which can be fit¬ 
ted with the Calabrese-Cardy formula, signaling the exis¬ 


tence of a critical (gapless) ground state, see Fig. SI 'b). 

Gapping mechanism. We now discuss the nature 
of the gapping mechanism responsible for the magnetic 
crystals studied in the text. The combined action of the 
Raman coupling 77 1 with 77int gives rise to two commut¬ 
ing Sine-Gordon terms: 

/ dx cos[x/2{q(j)c + &s)]+ / dx cos[V2{q^c-Ss)]. (S9) 



FIG. SI. Entanglement entropy S{£) of the ground state of a 
SU(2) fermionic gas for kso = rr/2 (a) and kso = tt/S (b); 
simulation parameters are v = 1/2, Q,/t = 1 and U/t —>■ oo. 
The fitted value of the central charge in panel (b) is c = 
1.05 ±0.15. 


Additionally, it is known that when z/ = 1, f, f,... the 
bare interaction 77int leads to the additional Mott terms 


TfMott.iy oc cos ( V 8 i/ V. 


c) , 


(SIO) 


which constitute a gapping mechanism for the charge de¬ 
grees of freedom when . We now discuss the 


interplay of the terms (S9) and (SIO) for different fillings 


v = l. 

When U = 0 and H 7 ^ 0 the ground state is fully 
gapped (band insulator) due to the presence of the 
terms ( |S9| ) with q = 1 (which are relevant for Kc < 3). If 
U > 0, the additional Mott term (SIOI with ix = 1 (which 
is relevant for Kc < 1) has to be taken into account. As 
the arguments of these three terms commute, (SIO) can¬ 


not modify the nature of the gapped phase or induce a 


phase transition (note also that the terms (S9) are more 


relevant). Furthermore, for U > 0 we have numerically 


checked that the gap induced by the terms (S9) is en¬ 
hanced with respect to the non-interacting regime [49j . 
as we can see by studying the correlation length ^ associ¬ 
ated to \{cl,m^3.Tn)\ ^ (see Fig. 

values of U. The gapped phase is stabi 


S2) for different 
ised by U since 
the correlation length ^ decreases (and thus the gap is 
enhanced) when U is increased. The additional fact that 
for fl = 0 and U ^ 0 the ground state is gapless (only 
charge degrees of freedom are gapped out) asserts that 
we are not observing a standard Mott insulator. 

13 = 1/3. When [7 = 0 and fl 7 ^ 0 the ground state 
is gapless. An incompressible phase is stabilised if H 7 ^ 0 
and V > 0 regardless the value of U. In this case we have: 


the terms (S9) with q = 3 (relevant if Kc < 1/3) and the 
Mott term (SIO) with = 1/3 (relevant if Kc < 1/9). 


Taking into account that Kc > 1/8, since we only have 
a nearest-neighbour interaction [7D], we obtain that the 
Mott term is not relevant and it cannot constitute a gap¬ 
ping mechanism. We thus conclude that the incompress¬ 
ible phase we observe is due to the interplay of interac¬ 
tions and magnetic field only. This is further supported 
by the observation that if we set H = 0 we can check 
numerically that the ground state is gapless and the en- 
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FIG. S2. Single-particle correlator \{c\,^Cj,m)\ for different 
values of the interaction strength U\ v = 1, fl/t = 1. 


tanglement entropy scales logarithmically with a central 
charge c = 2 (implying four gapless modes) for all the 
possible values of U and V. 


Supplementary Material for the case 1 = 1 


Numerical results for kso = ’’'/S and 'H 2 0. In ad¬ 

dition to the data for the fully-gapped phases for kso = 
tt/S discussed in the main text, in Fig. we present 
the density (fij) and magnetization profiles (M“) of the 
gapped phases for kso = 7 r /6 and fillings ly = 5 , 5 ,!. 


In Fig. S4 a) we also show the landscape of the on-site 

2tv\ 


eigenenergies 


= 217 cos 


27+1 


2 fcsoj 


(note that 


they explicitly depend on kso, and in the present case 
they display a six lattice site periodicity), in panels (b- 
d) we show the density profiles in the rotated spin basis. 





-(A7A=0> 


site(j) 


FIG. S4. Energy-spin structure Sj^\ for kso = ir/S (a). Den¬ 
sity plots (fij^x) for fillings n = ^ (b), | (c) and | (d) obtained 
through DMRG simulations (see the caption of Fig. S3 for the 
parameters). 



1/Vl 



FIG. S3. Gharge {fii} and magnetic order (M“) of the frac¬ 
tional phases at = I (a), | (b) and | (c) for kso = 7r/6 
as obtained from DMRG simulations of a system of length 
L = 96 with fl/t = 1. For the interaction parameters see the 
panels. Since the system is a crystal with small boundary ef¬ 
fects, we only plot the central part of the system for a better 
readability. 


FIG. S5. The oscillations of the currents (jfrn) vanish in the 
thermodynamic limit; a = (0.219 ± 0.001) and 6 = (5 • 10“® ± 
4- 10"®). 



i 

FIG. S6. DMRG simulation of the entanglement entropy of 
the state for I = 1, L — 192, i/ = j, D/i = 1, U/t —>■ 00 for 
kso = 7r/3. Thin red line is a fit with the Calabrese-Gardy 
formula which yields c = 1.0 ± 0.1. 
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Remarkably, they resemble a diluted version of the den¬ 
sity profiles presented in the main text, suggesting that, 
at a fixed u, almost the same physics is obtained by scal¬ 
ing kso and the number of fermions by the same factor. 
However, different (e.g. longer range) interactions may 
be necessary to stabilize phases with the same v. 

Numerical results for kso = ’’"/S and ‘H 2 = 0. We 
show that the oscillations of the chiral currents vanish 
in the thermodynamic limit, and can therefore be in¬ 
terpreted as a boundary/finite-size effect (see Fig. S5|. 
Conversely, the bulk value of such currents is indepen¬ 
dent from the system size. 

Gapless phases with 'H 2 7 ^ 0. Finally, we show that 
the range of the interactions is an essential ingredient 


to stabilise fully-gapped phases when 'H 2 7 ^ 0. Fig. ^ 
demonstrates that a simple contact interaction cannot 
stabilise a gapped phase at v = 1/4 with kso = tt/S: 
indeed, the ground-state von Neumann entanglement en¬ 
tropy displays a dependence on the subsystem size i 
which is typical of a gapless phase. 




